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The renin-angiotensin-aldosterone system (RAAS) plays a key role in the pathogenesis 
of cardiovascular disorders including hypertension and is one of the most important tar- 
gets for drugs. A whole body physiologically based pharmacokinetic (wb PBPK) model 
integrating this hormone circulation system and its inhibition can be used to explore the 
influence of drugs that interfere with this system, and thus to improve the understand- 
ing of interactions between drugs and the target system. In this study, we describe the 
development of a mechanistic RAAS model and exemplify drug action by a simulation of 
enalapril administration. Enalapril and its metabolite enalaprilat are potent inhibitors of the 
angiotensin-converting-enzyme (ACE). To this end, a coupled dynamic parent-metabolite 
PBPK model was developed and linked with the RAAS model that consists of seven coupled 
PBPK models for aldosterone, ACE, angiotensin 1, angiotensin 2, angiotensin 2 receptor 
type 1, renin, and prorenin.The results indicate that the model represents the interactions 
in the RAAS in response to the pharmacokinetics (PK) and pharmacodynamics (PD) of 
enalapril and enalaprilat in an accurate manner. The full set of RAAS-hormone profiles and 
interactions are consistently described at pre- and post-administration steady state as well 
as during their dynamic transition and show a good agreement with literature data. The 
model allows a simultaneous representation of the parent-metabolite conversion to the 
active form as well as the effect of the drug on the hormone levels, offering a detailed 
mechanistic insight into the hormone cascade and its inhibition. This model constitutes a 
first major step to establish a PBPK-PD-model including the PK and the mode of action 
(MoA) of a drug acting on a dynamic RAAS that can be further used to link to clinical 
endpoints such as blood pressure. 

Keywords: physiologically based pharmacokinetic model, cardiovascular, renin-angiotensin-aldosterone system, 
enalapril, enalaprilat 



INTRODUCTION 

Hypertension is a worldwide epidemic, which affects all ages 
and racial populations. It has a very high incidence and is the 
leading cause of cardiovascular mortality (Roger et al., 2012). 
There is still a lack of mechanistic models in the published 
literature that integrate the available physiological knowledge 
and relate the pharmacokinetics (PK) of hypertension drugs to 
their pharmacodynamics (PD) effects. One model addressing 
among others the RAAS and the blood pressure, is the Guy- 
ton's physiological model (Guyton et al, 1972), it was extended 
and published recently (Montani and Van Vliet, 2009; Osborn 
et al, 2009). Although this and other models may contain all 
the parameters listed in this present study and although they 
might have been validated against several physiological situa- 
tions they are not based on a generic whole body physiologi- 
cally based PK (wb PBPK) background model. Other approaches 



concerning the hormone cascade, as for example by White 
et al. (1989), Levitt and Schoemaker (2006), Guillaud and Han- 
naert (2010), Zhou et al. (2012), have different focuses or 
approaches. 

To overcome the lack of detailed physiologically based mod- 
els for cardiovascular PD, we aimed to establish a mechanistic 
model for the RAAS covering all relevant biological processes 
depicted in Figure 1. Renin is synthesized from the enzyme 
precursor prorenin, which is produced mainly in the kidney 
and secreted into the plasma by the granular cells (Krop et al, 
2008). The active enzyme renin is stored in the cells of the 
juxtaglomerular apparatus is released immediately on stimula- 
tion of these cells (Danser et al., 1998). After secretion into the 
plasma, renin converts the hepatically synthesized inactive hor- 
mone angiotensinogen (AGT) to angiotensin 1 (Angl). Angl is 
then converted in the plasma compartment by the membrane 
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bound angiotensin-converting-enzyme (ACE), expressed by the 
vascular endothelium, to angiotensin 2 (Ang2). Ang2 has an 
inhibiting effect on renin synthesis and secretion and thus on 
the plasma levels of Angl and Ang2 (Johns et al, 1990). Ang2 
is a vasoactive peptide that induces volume (and Na) retain- 
ing effects as well as vasoconstricting effects and an increase in 
blood volume, leading to increased blood pressure. It potenti- 
ates aldosterone secretion by interacting with its angiotensin 2 



receptor type 1 (ATI). Aldosterone increases the blood pressure 
as well (Gornall et al, 1960) by sodium and extracellular vol- 
ume retention. Several cardiovascular disease states are associated 
with changes in circulating Ang2 (Laragh et al., 1972; Laragh, 
1995). 

Inhibitors of the RAAS can be classified by their mode of action 

(Mo A; Cagnoni et al, 2010; Eisenberg, 2011; Aronson and Krum, 
2012) into five categories: 




FIGURE 1 | Schematic representation of the 
renin-angiotensin-aldosterone system (RAAS). Renin converts 
the hepatically secreted angiotensinogen (AGT) to angiotensin 1 
(Ang1). Ang1 is then converted to angiotensin 2 (Ang2) by the 
membrane bound angiotensin-converting-enzyme (ACE). Ang2 then 



binds to the angiotensin 2 receptor type 1 (ATD.The most important 
effects of AT1 binding are depicted either: the aldosterone secretion 
by the adrenals and the inhibition of renin activation from prorenin. 
Different processes are depicted by different line styles (see upper 
right box). 
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Table 1 | Parameters and variables used for the different actors in the RAAS. 



Parameter/variable name Value Unit Reference 





Molecular weight 


53.15 




kDa 




UniProt Consortium (2012) 




Plasma concentration* 


0.60 




jimol ■ L _1 




In agreement with Kobori et al. (2009) and Wang et al. 


(2012), 












consistent with the range given by Katsurada et al. (2007) 


fl/2 


11.29 




min 




Optimization 




^Secretion 


7.81 




L- min -1 




Optimization 




Molecular weight 


48 




kDa 




Sealey etal. (1980) 




Plasma concentration* 


3.62E- 


-07 


jimol ■ L _1 




Consistent with Juillerat et al. (1990), Nussberger et al 


(2002) 


^cat 


3.60 




min -1 




Optimization 






5.04E- 


-02 


jimol ■ L -1 




Optimization 




f 1/2 


1.44 




min 




Optimization 






Molecular weight 


180 




kDa 




Thevananther and Brecher (1999) 




k+ 


35.16 




L- ixmol -1 • 


min 1 


Results from IV Ena model 




k~ 


2.15E- 


-02 


min -1 




Results from IV Ena model 




kca\ 


4.48E- 


-04 


min -1 




Optimization 




K m 


4.12E- 


-03 


jimol ■ L -1 




Optimization 






Molecular weight 


1.30 




kDa 




Wang etal. (2012) 




Plasma concentration* 


7.91 E- 


-06 


jimol ■ L -1 




Consistent with Juillerat et al. (1990), Nussberger et al 


(2002) 


f 1/2 


0.72 




min 




Optimization 




Ang2 


Molecular weight 


1.05 




kDa 




Noda etal. (2006) 




Plasma concentration* 


4.84E- 


-06 


(j,mol ■ L -1 




Consistent with Juillerat et al. (1990), Nussberger et al 


(2002) 


Log P 


-1.7 








Wang etal. (2012) 




^ inhibition 


1.87E- 


-08 


(j,mol ■ L -1 




Optimization 




f 1/2 


1.54 




min 




Consistent with Boyd et al., 1967) 




Molecular weight 


41 




kDa 




UniProt Consortium (2012) 




Cp,ef 


7.69 




\kTT\0\ ■ L -1 




Optimization 




k + 


23.50 




L ■ |imol _1 ■ 


min 1 


Optimization 




k~ 


3.07E- 


-02 


min -1 




Optimization 






Molecular weight 


0.36 




kDa 




Wang etal. (2012) 




Plasma concentration* 


2.07E- 


-04 


\kTT\0\ ■ L -1 




Consistent with Juillerat et al. (1990), Nussberger et al 


(2002) 


Log P 


1.08 








Wishart etal. (2009) 




Protein binding 


50 




% 




Pardridge (2011) 




^Production 


1.51 E- 


-06 


(imol ■ min" 


-1 


Optimization 




^Secretion 


7.81 




L- min -1 




Optimization 




f 1/2 


12.87 




min 




Optimization 






Molecular weight 


57 




kDa 




Sealey et al. (1980) 




Plasma concentration* 


1.21 E- 


-06 


(j,mol ■ L -1 




Based onTu et al. (2012) with a concentration of 69 pg ■ 


ml- 1 


Renal plasma concentration intracellular 


1.1 E— ( 


34 


(j,mol ■ L -1 




Optimization 




^MaxKid 


1.62 




(imol ■ min" 


-1 . L -1 


Optimization 




Km Kid 


4.68E- 


-02 


limol ■ L -1 




Optimization 




^Secretion 


4.68E- 


-02 


L- min -1 




Optimization 




^Kidney (intracellular volume kidney) 


0.24 




L 




PK-Sim® 




f 1/2 


5.07 




min 




Optimization 






KmLiv 


710 




|xmol-L 1 




Tabata et al. (1990), Abu-Zahra and Pang (2000) 





'Plasma concentration in venous blood at pre-administration steady state. 
'Are assumed to circulate only in the plasma. 
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Table 2 | Physicochemical and PK properties of enalapril and enalaprilat. 





Enalapril 


Enalaprilat 


Unit 




Reference 




Molecular weight 


376.5* 


348.4* 


g ■ mol -1 




Wang etal. (2012) 




Log P 


0.07* 


-0.74* 


Log units 




Ranadive et al. (1992), Remko (2007) 


pKa 


3.74 (4.75*) 


2.03* 






Remko (2007) 








4.03* 










Plasma protein binding 


45 (50*) 


50-60 (50*) 


% 




Sirianni and Pang (1998), Knox et al. (2011) 


Hepatic clearance 


354 




ml ■ min -1 




Noormohamed et al. (1990) 


Renal clearance 


18 


8-9.5 


L-h" 1 




Dhareshwar (2007) 




* Values used in the model. 














Table 3 | Population specific information about subjects, optimized parameters, and model specific parameters. 






Biollazetal.(1982) 


Gu et al. (2004) 


Lee etal. (2003) 


Lu etal. (2009) Najib etal. (2003) 


Noormohamed 














et al. (1990) 


No. of individuals 


12 


20 


12 


20 


24 


12 


Male/female 


12/0 






20/0 


24/0 


12/0 


Age (years) 


22-33 








23.25±4.55 


21-25 (mean 22) 


Weight (kg) 










73.38 ±9.39 


56-80 (mean 70) 


Height (cm) 










175.96 ±766 




Healthy 


Yes 


Yes 


Yes 


Yes 


Yes 


Yes 


Dose (mg) 


10 


10 


20 


2-5* 


20 


20 


Fasted 


Yes 


No 


No 


Yes 


Yes 


Yes 




^maxLiv (|imol ■ min -1 -L _1 ) 


184.87 


178.18 


131.88 


148.05 


146.50 


194.36 


Renal clearance Ena (L- min -1 


■kg" 1 ) 6.02E-03 


4.51 E-3 


6.53E-03 


5.02E-03 


4.61 E-03 


5.50E-03 


Renal clearance Enaat (L- min" 


1 - kg" 1 ) 2.99E-04 


8.01 E-4 


4.99E-04 


6.30E-04 


6.37E-04 


6.21 E-04 


ACE ref cone. (CR e f) (|xmol ■ L~ 


') 2.59 


1.86 


4.32 


1.30 


2.14 


4.15 


PARAMETERS OBTAINED FROM PK-Sim® ACCORDING TO DIFFERENT AGE 










^Liver (U 


1.56 


1.57 


1.57 


1.35 


1.49 


1.42 



"Two enalapril maleate capsules were administered at the same time, containing 5mg enalapril each. 
' V mLl v values for carboxylesterase enzyme and the conversion from Ena to Enaat in Eg. 8. 



ACE inhibitors (ACEIs), 
Direct renin inhibitors (DRIs), 
Angiotensin-2-receptor-blockers (ARBs), 
Aldosterone synthase inhibitors (ASIs), and 
Potassium sparing diuretics. 

A common drug of the first type of inhibitors is enalaprilat (Enaat) , 
a potent and reversible ACE inhibitor (MacFadyen et al., 1993; 
Dhareshwar, 2007). Enaat is poorly absorbed from the gastroin- 
testinal tract. Therefore, an ethyl ester promoiety was incorporated 
within the administered parent drug enalapril (Ena) to increase 
its oral bioavailability. After absorption of inactive Ena by the 
peptide carrier-mediated transport system in the intestine, the 
prodrug Ena is converted to Enaat by the carboxylesterase fam- 
ily of enzymes activating it for ACE inhibition (MacFadyen et al, 
1993; Dhareshwar, 2007). 

MATERIALS AND METHODS 

A model of the RAAS based on available information was estab- 
lished (Eissing et al., 2011). To evaluate the RAAS model, a 



generic PBPK model of Ena and Enaat was developed and val- 
idated with clinical data that was then coupled to the RAAS 
model. 

SOFTWARE 

The RAAS model was created using the Systems Biology soft- 
ware suite that consists of the wb PBPK software tool PK-Sim® 
(Version 4.2, Bayer Technology Services, Leverkusen, Germany) 
and the general purpose modeling software MoBi® (Version 2.3, 
Bayer Technology Services, Leverkusen, Germany). This soft- 
ware has successfully been applied for diverse modeling stud- 
ies in laboratory animals and humans (Willmann et al., 2003, 
2004, 2007, 2010; Edginton et al, 2006a,b,c, 2008; Brochot et al, 
2010; Eissing et al., 2011). The parameter identifications for 
MoBi models were done using MATLAB® (R2010b, The Math- 
Works, Inc., Natick, MA, USA) and the MoBi Toolbox for MAT- 
LAB. Identification algorithm for Ena and Enaat models was 
Monte Carlo, while RAAS identification algorithm was frnin- 
search (fminsearch uses the Nelder-Mead Simplex method). The 
identification error function was in both cases error root mean 
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Table 4 | Reference expression values of ACE per organ as included in 
PK-Sim 5® expression database (Meyer et al., 2012) from Unigene 
(2010). 



Name 


ACE 


AT1 


Gene ID 


1636 


185 


Hs 


654434 


477887 


Organ 


Expression £ Ref 


in% (absolute 




value) of highest expression value 


Brain 


15 (12.96) 


5 (0.20) 


Fat 


6 (5.22) 


0(0) 


Gonads 


2 (2.01) 


4 (0.16) 


Heart 


0 (0.16) 


9 (0.35) 


Kidney 


0 (0.41) 


100 (3.85) 


Large Intestine 


2 (1.44) 


27 (1.04) 


Liver 


0(0) 


76 (2.94) 


Lung 


100 (86.38) 


6 (0.22) 


Muscle 


0 (0.12) 


8 (0.30) 


Pancreas 


0 (0.23) 


0(0) 


Plasma 


6 (5.55) 


0(0) 


Skin 


0(0) 


57 (2.19) 


Spleen 


0(0) 


89 (3.41) 


Small Intestine 


1 (0.32) 


27 (1.04) 


Stomach 


3 (2.86) 


0(0) 



square. For additional information, see the software manual or 
http://www.systems-biology.com/products. 

DEVELOPMENT OF THE RAAS-PBPK MODEL 

Physiological data presented in Table 1 contains the information 
gained from the literature or, if not available by optimization based 
on data for Angl, Ang2, aldosterone, renin, prorenin, ACE, ACT, 
and the ATI. Optimization was done with the overall coupled 
model and with all RAAS parameters (Table 1 ) simultaneously. For 
each of them an individual PBPK model was developed which con- 
tains many physiological background parameters (Eissing et al, 
201 1) that are derived from the few input parameters presented in 
Table 1. Mass transfer within the plasma (n) between organs is 
calculated as organ specific blood flow ( Q or gan) times the drug con- 
centration within the plasma of the respective organ (C p i s _ orga n) 
times the overall plasma fraction (1 — HCT) which is one sub- 
tracted by the hematocrit (HCT) value. Passive transport into the 
interstitial and intracellular space is realized by multiplying organ 
specific partition coefficients K m% . 



1 — Qorgan ' ( 1 H CT) • Cpl s _ or g an 



(1) 



The interactions in the coupled model were described as 
depicted in Figure 1 reflecting the biological processes as detailed 
in the following. ACT, prorenin, renin, and aldosterone are the 
four hormones and enzymes synthesized in cells and then secreted 
endocrinally. It was assumed that these hormones and enzymes 
are produced intracellular and secreted into the plasma by a 
first order rate (r 2 in |xmol-min _1 ) proportional to the intra- 
cellular concentration (Ccdiular Hormone) m the respective organ 



with the rate constant /cs eC retion (in L-min -1 ). The renin secre- 
tion r$ (|xmol • min -1 ) is described by a "competitive" inhibition 
equation including the conversion of prorenin to renin influ- 
enced by Ang2 (see Figure 1), describing the conversion from 
prorenin to renin, and simultaneously describing the secretion 
from intracellular to plasma compartment of the kidney: 



''3 



VkaxKid ' Vkidney ' Cp ri 



(0 



Cp r , 



(t) + KrnKid 



/ 1 + /Q32ffl\\ 

\ \ ^inhibition / / 



(2) 



^mKid an d ^MaxKid are the respective enzymatic activity con- 
stants, Vjddney is me intracellular volume of the kidney, Cp r0 renin is 
the renal intracellular prorenin concentration (parameters under 
prorenin in Table 1). inhibition represents the inhibitory constant 
for Ang2 and CAng2 is the renal plasma concentration of Ang2. The 
degradation rate r$ ([xmol • min -1 • L _1 ) is described by Guillaud 
and Hannaert (2010): 



U = C (f) 



ln2\ 
J1/2) 



(3) 



fi/2 determines the half-life and C(t) being the respective plasma 
concentration over time. These half-lives were optimized for the 
hormones and enzymes (see Table 1). The conversion of ACT 
to Angl and from Angl to Ang2 take place in the plasma com- 
partments of all organs (Eissing et al, 2011) and is an enzymati- 
cally regulated process described by a Michaelis-Menten equation 
(Gould et al, 1980; Ehlers and Kirsch, 1988). The change of hor- 
mone mass over time (rs in |xmol • min -1 ) then depends on the 
enzyme concentration (Co) and the prohormone (Cp ro hormone) 
concentration in the respective organs, as well as on the enzyme 
activity represented by the K m and fc cat of the enzyme (Meyer et al, 
2012). This can be expressed by: 



Vo • Co (t) ■ Cprohormone (0 
^Prohormone (0 + f^m 



(4) 



The calculation refers to the respective organ plasma volumes 
Vq in L (taken from PK-Sim® database). ACE is most notably 
bound to the endothelium (Dzau et al, 2001) and thus not circu- 
lating in the model. The concentration of ACE and ATI in each 
organ (Co) is calculated via a reference concentration (CR e f, see 
Table 3) times the percentage of reference expression (ER e f) from 
literature (see Table 4, factor ER e f =100 for the highest expression 
value): 



Q = CR e f 



^Ref 

100 



(5) 



Ang2 binds to the ATI receptor, which influences the aldosterone 
secretion. ATI binding rate ( r-j in (imol • min~ 1 • IT 1 ) is described 
by a typical equation with k + and k~ binding constants for asso- 
ciation and dissociation and the respective concentrations for ATI 
(C 0 ), Ang2 (C 2 ), and the AT1-Ang2 complex (C Co mplex): 



r 7 = (k+ ■ Co • (f) • C 2 (t)) - (k_ ■ C comp i ex (0) 



(6) 
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Time [min] Simulated (Enaat) cone. [uM] 



• 


Adolescent data 


Simulation 


▼ 


Elderly data 


Regression line 



FIGURE 2 | (A) Plasma concentration-time profiles of Enaat after IV observed plot for the same simulation with the Hockings data are 

Enaat administration. The solid blue line is the simulation result shown, indicated as points (adolescents) and triangles (elderly). The solid 

compared to the data of Hockings et al. (1986) which are indicated as black regression line denotes the line of identity (observed outcome is 

points (adolescents) and triangles (elderly). (B)The simulated versus equal to simulated outcome). 



The ATI concentrations were calculated as explained in Eq. 
6, the respective relative expression values (£Ref) are shown in 
Table 4, and the reference concentration (CR e f) is presented in 
Table 1. As shown in Figure 1 ATI is membrane bound too and 
thereby does not circulate in the plasma. Aldosterone synthesis ( r& 
in ixmol- min -1 ) in our model was directly linked to the com- 
plex formation of Ang2 bound to the ATI ( Complex )> which is 
modeled as 

r 8 — (^secretion * C CO mpl e x (0) ~t~ ^production (7) 



concentrations after IV Enaat (Hockings et al., 1986) administra- 
tion. The coupling of parent and metabolite was realized with a 
hydrolytic reaction in the hepatocytes (Dhareshwar, 2007). This 
reaction is mediated by different carboxylesterase enzymes. K m 
values for metabolism to Enaat in Table 1 were taken from mea- 
surements in rat hepatocytes (Tabata et al, 1990; Abu-Zahra and 
Pang, 2000). The conversion rate {rg in |xmol • min -1 ) from Ena 
into Enaat in the hepatocytes with the respective intracellular 
volume Vu vel - (see Table 3), the concentrations of Ena (CE na ) 
and Enaat (CEnaat) is modeled by Michaelis-Menten kinetics: 



The secretion is thereby described to be dependent on intracellu- 
lar concentration of Ang2 receptor complex multiplied by the rate 
constant ^secretion (in L-min -1 ) and additionally represents the 
continuous secretion as fcp ro d uct i on (Veldhuis et al, 2008). 

DEVELOPMENT OF THE ENALAPRIL/ENALAPRILAT PBPK MODEL 

A PBPK model for oral Ena including the metabolism to active 
Enaat was developed first and used as test case for the dynamic 
behavior of the RAAS model. Based on intravenous (IV) data, 
parameters for Enaat were determined and then used in the parent- 
metabolite model. The required physicochemical input properties 
and clearance processes for the model are presented in Table 2. 
The model simulations for oral Ena were compared to differ- 
ent experimental data (Biollaz et al, 1982; Noormohamed et al., 
1990; Nussberger et al., 2002; Lee et al, 2003; Najib et al, 2003; 
Gu et al, 2004; Lu et al, 2009). The six different subject and 
dataset properties (all included to the respective simulations) are 
shown in Table 3. Only the five parameters presented in Table 3 
were optimized with MATLAB®. The respective IV Enaat sim- 
ulations were compared to experimental data of Enaat plasma 



Liv ' Vuver • CEna (0 
rg = (8) 

CEna (f) + ^mLiv 

The respective parameters for Michaelis-Menten kinetics are 
VMaxLiv and K m u v . 

ASSEMBLING OF THE MODELS TO A FULL DYNAMIC PBPK-RAAS 
MODEL 

To couple the drug model to the full RAAS model the MoA 
was included. Enaat binds to ACE in the RAAS (see Table 1, 
ACE parameters) generating an inactive complex. The forma- 
tion rate of Enaat- ACE complex (Ccomplex) is realized as pre- 
sented in equation 6 for all plasma compartments. ACE con- 
centration is based on the relative expression value in Table 4 
which was obtained from Unigene (2010) and processed as 
described in Meyer et al. (2012). The inhibition of ACE occurs 
by ACE binding to Enaat, forming an ACE-Enaat- complex, 
thereby blocking ACE for the interaction with Angl. The 
binding is implemented in the plasma compartment of each 
organ. 
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Time [min] 

B 




10" 3 10~ 2 10" 1 10 1 



Simulated (Ena) cone. [uM] 



A) 




Ena simulation 




• 


Ena data 




— 


Enaat simulation 




▼ 


Enaat data 


B) 


• 


Biollaz et al. 


C) 




Gu et al. 






Lee et al. 




X 


Lu et al. 




X 


Najib et al. 




+ 


Noormohamed et al. 




Simulated (Enaat) cone. [uM] 



FIGURE 3 | (A) Exemplary plasma concentration-time curves of Ena and 
Enaat from the coupled Ena-Enaat model after oral Ena application compared 
to experimental data (Gu et al., 2004). The red solid line is the simulated Ena 
plasma concentration, while the blue solid line is the simulation result of 
Enaat. The respective data of Ena are indicated as circles, of Enaat as 
triangles. (B) Simulated versus observed plot for Enalapril from the coupled 
model, compared to the data by Biollaz et al. (1982), Noormohamed et al. 



(1990), Lee et al. (2003), Najib et al. (2003), Gu et al. (2004), Lu et al. (2009), 
with SD. The solid black regression line denotes the line of identity (observed 
outcome is equal to simulated outcome). (C) Predicted versus observed plot 
for Enalaprilat from the coupled Ena-Enaat model, compared to the data by 
Biollaz et al. (1982), Noormohamed et al. (1990), Lee et al. (2003), Najib et al. 
(2003), Gu etal. (2004), Lu et al. (2009). The solid black regression line denotes 
the line of identity (observed outcome is equal to simulated outcome). 



RESULTS 

ORAL ADMINISTRATION OF ENALAPRIL 

A PBPK model of IV Enaat has been developed with the data 
indicated in Table 2. In Figure 2 the resulting PK-simulations for 
IV Enaat are shown and the goodness-of-fit of the model based 
on data is depicted (Hockings et al, 1986). This model was used 
to determine the parameters related to Enaat, which are log P, 
Molecular weight, pKa, plasma protein binding and clearance (see 
Table 2). Additionally the constants of ACE binding were iden- 
tified (see Table 1). The constant for association is much higher 
than the dissociation constant, as stated by Reynolds (1984). This 
parameterization was used within the coupled model without 
changes. 

Resulting plasma concentrations of ACE based on the explained 
calculation (peripheral venous plasma concentration in the model 
validated by data from Gu et al. (2004): 7.17E-4 umol ■ L _1 ) 
correspond nicely with literature data of 2.52E— 3 (imol- L _1 by 
Brice et al. (1995). After integrating the metabolite model into the 
oral model, the coupled model exhibits the behavior presented 
in Figure 3. Including the population dependent parameters for 



the different publications (see Table 3), the six model simulations 
show excellent agreement with the data. The serum concentration 
profile of Enaat exhibits a prolonged terminal phase, apparently 
representing the small fraction of the administered dose that has 
been bound to ACE. As obvious from Figure 3, the model Ena 
simulations and Enaat predictions accurately fit the observed data. 

RAAS INHIBITION BY ENALAPRILAT 

The PBPK model for the circulating, endocrine RAAS was devel- 
oped to describe the steady state described in literature. It corre- 
sponds to the pre-administration steady state discussed further 
in the following. The model for Ena was then coupled to the 
RAAS model to explore the dynamical behavior of the system 
under inhibition. In Figure 4 the results for the full RAAS model 
under administration of 10 mg Ena are shown at t = 1500 min 
(administration at t = 1500 indicated by the red arrow). All para- 
meters of the Ena-Enaat model were used without changes. The 
pre- and post-administration steady states as well as the dynam- 
ical transition after administration at t = 1500 min (red arrow) 
show excellent agreement to literature data (Table 1, Figure 4A). 
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— Aldost simulation 


— Ang2 simulation 


• Aldost data 


■ Ang2 data 


— Ang1 simulation 


— Renin simulation 


▼ Ang1 data 


x Renin data 


— Regression line 


— ► Timepoint of administration 



FIGURE 4 | (A) Plasma concentration-time curves of hormones 
included into the RAAS. Hormones are in a steady state until 
administration of oral Ena at t = 1500 min indicated by the red arrow. 
The respective data are indicated by different symbols (Nussberger 
et al., 2002). (B) Simulated versus observed plasma concentrations 



of different actors of the RAAS before and after inhibition of ACE by 
Ena. The different hormones are depicted by different symbols 
(Nussberger et al., 2002) with their SD. The solid black regression 
line denotes the line of identity (observed outcome is equal to 
simulated outcome). 



In detail, Figure 4 shows the MoA of Ena. After the administra- 
tion of Ena, ACE is inhibited. As can be seen, the conversion of 
Angl to Ang2 is inhibited, lowering the Ang2 and aldosterone 
plasma levels. The negative feedback loop from the Ang2 to the 
renin synthesis is inhibited, leading to an increase of the renin 
plasma concentration, and thereby the Angl levels. The compar- 
ison of these induced changes in plasma levels by a scatter plot 
(Figure 4B) to study data from Nussberger et al. (2002) shows 
good agreement and corresponds also to older experimental data 
(Given et al, 1984; Johnston et al, 1984; Juillerat et al, 1990; Azizi 
etal., 1997). 

The simulation of exemplary Ena administration shows, 
that the full dynamical behavior of inhibition and the other 
RAAS-hormone levels can be represented accurately (see 
Figure 4). 

DISCUSSION 

A full mechanistic model of the circulating RAAS was developed. 
The model represents the steady state and dynamic behavior of 
the RAAS in an accurate manner. This was demonstrated by sim- 
ulating the effects of representative drug acting on the system. The 
Angl concentration is influenced by different reactions, Angl gen- 
eration, Ang2 generation, degradation of Angl, which lead to the 
pre-administration steady state. So it is rather difficult to evalu- 
ate the condition of substrate excess for the second, since the new 
generation of Angl is fast. In addition the reference concentra- 
tions Co are artificial constructs in a product with k cat , where both 
and K m are fitted parameters. Thus Michaelis-Menten equations 
were used, although the Co is much higher than the Angl con- 
centration for the steady state, as it was used for these enzymes 



to catalyze this reaction in the described way (Ehlers and Kirsch, 
1988; Guillaud and Hannaert, 2010). Although we find an excellent 
correlation of plasma levels and a good correlation of parameters 
to the articles cited before, it has not been proven in this arti- 
cle, that the optimized parameter values are unique to reach the 
reference data, since not for all parameters measured values are 
available. The Ena model was validated independently from the 
parameterization of the RAAS and can thus also serve as a valida- 
tion for the RAAS. Following the development of the Ena model 
based on literature data, subsequently only population dependent 
parameters were changed (listed in Table 3) to represent physiolog- 
ical differences in the studies. After coupling to the RAAS model 
the overall model was able to represent detailed information of 
all main components and predict circulating plasma concentra- 
tions of all included hormones by one coupled PBPK/PD model. 
The excellent representation of the hormone levels at steady state 
and during inhibition, implicate a good representation of the 
underlying physiological processes in the RAAS. 

The detailed physiological representation of relevant processes 
spanning from wb distribution to molecular interactions allows 
for, e.g., a mechanistic implementation of independent knowl- 
edge on drug-target interactions, and thereby an integration and 
translation of available knowledge. 

Thus the model constitutes a first major step to establish a 
PBPK-PD-model for drugs acting on blood pressure and heart 
rate. Generally blood pressure and heart rate represent the state 
of the art endpoint in clinical studies. Currently blood pressure 
is not yet included in the model. The inclusion of inhibitors 
from other classes of MoA's can be used to explore pharmaco- 
dynamic effects in the same model framework. Further extensions 
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of this model and an inclusion of the dynamic behavior of the 
blood pressure in response to changes in RAAS would allow 
a full mechanistic representation of the PBPK-PD relationship 
in cardiovascular diseases. Such a model could thus rational- 
ize treatment decisions in the area of cardiovascular diseases. 



By inclusion of a full process description from administration 
of a drug to its effect, i.e., the lowering of blood pressure, 
such a model allows to support cardiovascular drug research 
from target identification and validation to dosing regimen 
decision. 



REFERENCES 

Abu-Zahra, T. N., and Pang, K. S. 
(2000). Effect of zonal transport 
and metabolism on hepatic removal: 
enalapril hydrolysis in zonal, isolated 
rat hepatocytes in vitro and cor- 
relation with perfusion data. Drug 
Metab. Dispos. 28, 807-813. 

Aronson, D., and Krum, H. (2012). 
Novel therapies in acute and chronic 
heart failure. Pharmacol. Ther. 23, 
23. 

Azizi, M., Guyene, T. T., Chatellier, 
G., Wargon, M., and Menard, J. 
(1997). Additive effects of losartan 
and enalapril on blood pressure and 
plasma active renin. Hypertension 
29, 634-640. 

Biollaz, J., Schelling, J. L., Jacot 
Des Combes, B., Brunner, D. B., 
Desponds, G., Brunner, H. R., 
et al. (1982). Enalapril maleate 
and a lysine analogue (MK-521) 
in normal volunteers; relation- 
ship between plasma drug lev- 
els and the renin angiotensin sys- 
tem. Br. }. Clin. Pharmacol. 14, 
363-368. 

Boyd, G. W., Landon, J., and Peart, 
W. S. (1967). Radioimmunoassay 
for determining plasms-levels of 
angiotensin II in man. Lancet 2, 
1002-1005. 

Brice, E. A., Friedlander, W., Bateman, E. 
D., and Kirsch, R. E. (1995). Serum 
angiotensin-converting enzyme 
activity, concentration, and specific 
activity in granulomatous intersti- 
tial lung disease, tuberculosis, and 
COPD. Chest 107, 706-710. 

Brochot, A., Zamacona, M., and 
Stockis, A. (2010). Physio- 
logically based pharmaco- 
kinetic/ pharmacodynamic 
animal-to-man prediction of 
therapeutic dose in a model of 
epilepsy. Basic Clin. Pharmacol. 
Toxicol. 106, 256-262. 

Cagnoni, E, Njwe, C. A., Zaninelli, A., 
Ricci, A. R., Daffra, D., D'Ospina, A., 
et al. (2010). Blocking the RAAS at 
different levels: an update on the use 
of the direct renin inhibitors alone 
and in combination. Vase. Health 
Risk Manag. 6, 549-559. 

Danser, A. H.,Derkx, F. H., Schalekamp, 
M. A., Hense, H. W., Riegger, G. A., 
and Schunkert, H. (1998). Determi- 
nants of interindividual variation of 
renin and prorenin concentrations: 
evidence for a sexual dimorphism 



of (pro)renin levels in humans. /. 
Hypertens. 16, 853-862. 

Dhareshwar, S. S. (2007). "Case study: 
enalapril: a prodrug of enalaprilat," 
in Prodrugs, eds. V. J. Stella, R. T. Bor- 
chardt, M. J. Hageman, R. Oliyai, H. 
Maag, and J. W. Tilley (New York: 
Springer), 1221-1229. 

Dzau, V. J., Bernstein, K., Celermajer, D., 
Cohen, J., Dahlof, B., Deanfield, J., 
et al. (2001). The relevance of tis- 
sue angiotensin-converting enzyme: 
manifestations in mechanistic and 
endpoint data. Am. J. Cardiol 88, 
1L-20L. 

Edginton, A. N., Schmitt, W., Voith, B., 
and Willmann, S. (2006a). A mech- 
anistic approach for the scaling of 
clearance in children. Clin. Pharma- 
cokinet. 45, 683-704. 

Edginton, A. N., Schmitt, W., and Will- 
mann, S. (2006b). Application of 
physiology-based pharmacokinetic 
and pharmacodynamic modeling 
to individualized target-controlled 
propofol infusions. Adv. Ther. 23, 
143-158. 

Edginton, A. N., Schmitt, W., and Will- 
mann, S. (2006c). Development and 
evaluation of a generic physiologi- 
cally based pharmacokinetic model 
for children. Clin. Pharmacokinet. 
45, 1013-1034. 

Edginton, A. N., Theil, F. P., Schmitt, 
W., and Willmann, S. (2008). Whole 
body physiologically-based pharma- 
cokinetic models: their use in clini- 
cal drug development. Expert Opin. 
Drug Metab. Toxicol. 4, 1143-1152. 

Ehlers, M. R., and Kirsch, R. E. (1988). 
Catalysis of angiotensin I hydrolysis 
by human angiotensin-converting 
enzyme: effect of chloride and pH. 
Biochemistry 27, 5538-5544. 

Eisenberg, J. M. (2011). ACEIs, ARBs, 
or DR1 for Adults with Hypertension. 
Rockville, MD: Agency for Health- 
care Research and Quality. 

Eissing, T., Kuepfer, L., Becker, C, Block, 
M., Coboeken, K., Gaub, T., et al. 
(2011). A computational systems 
biology software platform for multi- 
scale modeling and simulation: inte- 
grating whole-body physiology, dis- 
ease biology, and molecular reac- 
tion networks. Front. Physiol. 2:4. 
doi: 10.3389/fphys.20 1 1 .00004 

Given, B. D., Taylor, T, Hollenberg, 
N. K., and Williams, G. H. (1984). 
Duration of action and short-term 
hormonal responses to enalapril 



(MK 421) in normal subjects. /. 
Cardiovasc. Pharmacol. 6, 436-441. 

Gornall, A. G., Grundy, H. M., and 
Koladich, C. J. (1960). Studies on 
aldosterone. III. Chronic effect on 
the blood pressure of rats. Can. }. 
Biochem. Physiol. 38, 43-50. 

Gould, A. B., Dewolf, R., Goodman, 
S., Onesti, G., and Swartz, C. 
(1980). Kinetic studies of the human 
renin and human substrate reaction. 
Biochem. Med. 24, 321-326. 

Gu, Q., Chen, X., Zhong, D., and Wang, 
Y. (2004). Simultaneous determi- 
nation of enalapril and enalapri- 
lat in human plasma by liq- 
uid chromatography-tandem mass 
spectrometry. /. Chromatogr. BAna- 
lyt. Technol. Biomed. Life Sci. 813, 
337-342. 

Guillaud, E, and Hannaert, P. (2010). A 
computational model of the circu- 
lating renin-angiotensin system and 
blood pressure regulation. Acta Bio- 
theor. 58, 143-170. 

Guyton, A. C, Coleman, T. G., and 
Granger, H. J. (1972). Circulation: 
overall regulation. Annu. Rev. Phys- 
iol. 34, 13-46. 

Hockings, N., Ajayi, A. A., and Reid, 
J. L. (1986). Age and the pharma- 
cokinetics of angiotensin convert- 
ing enzyme inhibitors enalapril and 
enalaprilat. Br. J. Clin. Pharmacol. 
21,341-348. 

Johns, D. W., Peach, M. J., Gomez, 
R. A., Inagami, T, and Carey, R. 
M. (1990). Angiotensin II regulates 
renin gene expression. Am. J. Physiol. 
259, F882-F887. 

Johnston, C. I., Jackson, B. J., Lar- 
mour, I., Cubella, R., and Casley, 
D. (1984). Plasma enalapril levels 
and hormonal effects after short- 
and long-term administration in 
essential hypertension. Br. J. Clin. 
Pharmacol. 18(Suppl.2),233S-239S, 
24 IS. 

Juillerat, L., Nussberger, J., Menard, 
J., Mooser, V., Christen, Y, Wae- 
ber, B., et al. (1990). Determinants 
of angiotensin II generation dur- 
ing converting enzyme inhibition. 
Hypertension 16, 564-572. 

Katsurada, A., Hagiwara, Y, Miyashita, 
K., Satou, R., Miyata, K., Ohashi, N., 
et al. (2007). Novel sandwich ELISA 
for human angiotensinogen. Am. J. 
Physiol. Renal Physiol. 293, 6. 

Knox, C, Law, V., Jewison, T, Liu, 
P., Ly, S., Frolkis, A., et al. (2011). 



DrugBank 3.0: a comprehensive 
resource for "omics" research on 
drugs. Nucleic Acids Res. 39, 8. 

Kobori, H., Alper, A. B. Jr., Shenava, R., 
Katsurada, A., Saito, T, Ohashi, N., 
et al. (2009). Urinary angiotensino- 
gen as a novel biomarker of the 
intrarenal renin- angiotensin sys- 
tem status in hypertensive patients. 
Hypertension 53, 344-350. 

Krop, M., De Bruyn, J. H., Derkx, 
F. EL, and Danser, A. H. (2008). 
Renin and prorenin disappearance 
in humans post-nephrectomy: evi- 
dence for binding? Front. Biosci. 13, 
3931-3939. 

Laragh, J. H. (1995). Renin- 
angiotensin-aldosterone system 
for blood pressure and electrolyte 
homeostasis and its involvement in 
hypertension, in congestive heart 
failure and in associated cardiovas- 
cular damage (myocardial infarction 
and stroke). /. Hum. Hypertens. 9, 
385-390. 

Laragh, J. EL, Baer, L., Brunner, H. 
R., Buhler, E R., Sealey, J. E., and 
Vaughan, E. D. Jr. (1972). Renin, 
angiotensin and aldosterone system 
in pathogenesis and management of 
hypertensive vascular disease. Am. J. 
Med. 52,633-652. 

Lee, J., Son, J., Lee, M., Lee, K. T, 
and Kim, D. H. (2003). Simulta- 
neous quantitation of enalapril and 
enalaprilat in human plasma by 96- 
well solid-phase extraction and liq- 
uid chromatogr aphy/tandem mass 
spectrometry. Rapid Commun. Mass 
Spectrom. 17, 1157-1162. 

Levitt, D. G., and Schoemaker, R. 
C. (2006). Human physiologi- 
cally based pharmacokinetic model 
for ACE inhibitors: ramipril and 
ramiprilat. BMC Clin. Pharmacol. 
6:1. doi:10.1186/1472-6904-6-l 

Lu, S., Jiang, K., Qin, E, Lu, X., and Li, 
F. (2009). Simultaneous quantifica- 
tion of enalapril and enalaprilat in 
human plasma by high-performance 
liquid chromatography-tandem 
mass spectrometry and its applica- 
tion in a pharmacokinetic study. /. 
Pharm. Biomed. Anal. 49, 163-167. 

MacFadyen, R. J., Meredith, P. 
A., and Elliott, H. L. (1993). 
Enalapril clinical pharmacoki- 
netics and pharmacokinetic - 
pharmacodynamic relationships. 
An overview. Clin. Pharmacokinet. 
25,274-282. 



w w w.f ro n t i ersin.org 



February 2013 | Volume 4 | Article 4 | 9 



Claassen et al. 



PBPK model of the RAAS 



Meyer, M., Schneckener, S., Ludewig, 

B. , Kuepfer, L., and Lippert, J. 
(2012). Using expression data for 
quantification of active processes in 
physiologically-based pharmacoki- 
netic modeling. DrugMetab. Dispos. 
40,892-901. 

Montani, J. P., and Van Vliet, B. N. 
(2009). Understanding the contri- 
bution of Guyton's large circula- 
tory model to long-term control of 
arterial pressure. Exp. Physiol. 94, 
382-388. 

Najib, N. M., Idkaidek, N., Adel, A., 
Admour, I., Astigarraga, R. E., De 
Nucci, G., et al. (2003). Bioequiv- 
alence evaluation of two brands of 
enalapril 20 mg tablets (Narapril and 
Renitec) in healthy human volun- 
teers. Biopharm. Drug Dispos. 24, 
315-320. 

Noda, T., Yaginuma, T., O'Rourke, M. 
P., and Hosoda, S. (2006). Effects 
of nifedipine on systemic and pul- 
monary vascular impedance in sub- 
jects undergoing cardiac catheteriza- 
tion. Hypertens. Res. 29, 505-513. 

Noormohamed, F. H., McNabb, W. R., 
and Lant, A. F. (1990). Pharmacoki- 
netic and pharmacodynamic actions 
of enalapril in humans: effect of 
probenecid pretreatment. /. Phar- 
macol. Exp. Ther. 253, 362-368. 

Nussberger, J., Wuerzner, G, Jensen, 

C, and Brunner, H. R. (2002). 
Angiotensin II suppression in 
humans by the orally active renin 
inhibitor Aliskiren (SPP100): com- 
parison with enalapril. Hypertension 
39.E1-E8. 

Osborn, J. W., Averina, V. A., and Fink, 
G. D. (2009). Current computational 
models do not reveal the impor- 
tance of the nervous system in long- 
term control of arterial pressure. 
Exp. Physiol. 94, 389-396. 

Pardridge, W. M. (2011). Targeted deliv- 
ery of hormones to tissues by 
plasma proteins. Compr. Physiol. 
335-382. 

Ranadive, S. A., Chen, A. X., and 
Serajuddin, A. T. (1992). Relative 



lipophilicities and structural- 
pharmacological considerations 
of various angiotensin-converting 
enzyme (ACE) inhibitors. Pharm. 
Res. 9, 1480-1486. 

Remko, M. (2007). Acidity, lipophilic- 
ity, solubility, absorption, and polar 
surface area of some ACE inhibitors. 
Chem. Pap. 61, 133-141. 

Reynolds, C. H. (1984). Kinetics of 
inhibition of angiotensin convert- 
ing enzyme by captopril and by 
enalapril diacid. Biochem. Pharma- 
col. 33, 1273-1276. 

Roger, V. L., Go, A. S., Lloyd-Jones, 
D. M., Benjamin, E. J., Berry, J. D., 
Borden, W. B., et al. (2012). Heart 
disease and stroke statistics - 2012 
update: a report from the American 
Heart Association. Circulation 125, 
e2-e220. 

Sealey, J. E., Atlas, S. A., and Laragh, J. 
H. (1980). Prorenin and other large 
molecular weight forms of renin. 
Endocr. Rev. 1, 365-391. 

Sirianni, G. L., and Pang, K. S. 

(1998) . Intracellular and not intra- 
luminal esterolysis of enalapril 
in kidney. Studies with the sin- 
gle pass perfused nonfiltering rat 
kidney. Drug Metab. Dispos. 26, 
324-331. 

Tabata, S., Yamazaki, H., Ohtake, Y, 
and Hayashi, S. (1990). Rapid 
conversion of the new angiotensin 
converting enzyme inhibitor 
ramipril to its active metabolite 
in rats. Arzneimittelforschung 40, 
865-867. 

Thevananther, S., and Brecher, A. S. 

( 1 999) . Isolation of angiotensin con- 
verting enzyme (ACE) binding pro- 
tein from human serum with an 
ACE affinity column. Can. J. Physiol. 
Pharmacol. 77, 216-223. 

Tu, W., Eckert, G. J., Pratt, J. H., and 
JanDanser.A. H. (2012). Plasma lev- 
els of prorenin and renin in blacks 
and whites: their relative abundance 
and associations with plasma aldos- 
terone concentration. Am. /. Hyper- 
tens. 25, 1030-1034. 



Unigene, (2010). National Center 
for Biotechnology Informa- 
tion (NCBI). Available at: 
www.ncbi.nlm.nih.gov/unigene 

UniProt Consortium, S. (2012). Reorga- 
nizing the protein space at the Uni- 
versal Protein Resource (UniProt). 
Nucleic Acids Res. 40, 18. 

Veldhuis, J. D., Keenan, D. M., and 
Pincus, S. M. (2008). Motivations 
and methods for analyzing pulsatile 
hormone secretion. Endocr. Rev. 29, 
823-864. 

Wang, Y, Xiao, J., Suzek, T. O., Zhang, J., 
Wang, J., Zhou, Z., et al. (2012). Pub- 
Chem's BioAssay Database. Nucleic 
Acids Res. 40, 2. 

White, K. P. Jr., Radke-Sharpe, N. 
E, Kaiser, D. L., and Owens, G. 
K. (1989). A dynamic model of 
angiotensin II infusion exper- 
iments. /. Biomed. Eng. 11, 
63-71. 

Willmann, S., Hohn, K., Edginton, A., 
Sevestre, M., Solodenko, J., Weiss, 
W., et al. (2007). Development of a 
physiology-based whole-body pop- 
ulation model for assessing the influ- 
ence of individual variability on 
the pharmacokinetics of drugs. /. 
Pharmacokinet. Pharmacodyn. 34, 
401-431. 

Willmann, S., Schmitt, W., Keldenich, 
J., and Dressman, J. B. (2003). 
A physiologic model for simulat- 
ing gastrointestinal flow and drug 
absorption in rats. Pharm. Res. 20, 
1766-1771. 

Willmann, S., Schmitt, W., Keldenich, 
J., Lippert, J., and Dressman, J. B. 
(2004). A physiological model for 
the estimation of the fraction dose 
absorbed in humans. /. Med. Chem. 
47, 4022-4031. 

Willmann, S., Thelen, K., Becker, C, 
Dressman, J. B., and Lippert, J. 
(2010). Mechanism-based predic- 
tion of particle size-dependent 
dissolution and absorption: cilosta- 
zol pharmacokinetics in dogs. 
Eur. J. Pharm. Biopharm. 76, 
83-94. 



Wishart, D. S., Knox, C, Guo, A. 
C, Eisner, R., Young, N., Gau- 
tam, B., et al. (2009). HMDB: 
a knowledgebase for the human 
metabolome. Nucleic Acids Res. 37, 
25. 

Zhou, X., Shang, D., Zhang, T., 
Li, L., Zhou, T., and Lu, W. 
(2012). Modeling of angiotensin 
II-angiotensin-(l-7) counter- 
balance in disease progression 
in spontaneously hyperten- 
sive rats treated with/without 
perindopril. Pharmacol. Res. 
9, 9. 

Conflict of Interest Statement: Karina 
Claassen is funded by and Thomas Eiss- 
ing, Michael Block, and Stefan Will- 
mann are employes of Bayer Tech- 
nology Services GmbH, the company 
that owns and commercializes the soft- 
ware platform used for the simulations 
described in the manuscript (PK-Sim® 
and MoBi®), as well as parent company 
stock owners. 

Received: 29 July 2012; accepted: 08 Jan- 
uary 201 3; published online: 08 February 
2013. 

Citation: Claassen K, Willmann S, 
Eissing T, Preusser T and Block 
M (2013) A detailed physiologically 
based model to simulate the phar- 
macokinetics and hormonal pharma- 
codynamics of enalapril on the cir- 
culating endocrine renin-angiotensin- 
aldosterone system. Front. Physio. 4:4. 
doi: 10.3389/fphys.2013.00004 
This article was submitted to Frontiers in 
Computational Physiology and Medicine, 
a specialty of Frontiers in Physiology. 
Copyright © 2013 Claassen, Willmann, 
Eissing, Preusser and Block. This is an 
open-access article distributed under the 
terms of the Creative Commons Attribu- 
tion License, which permits use, distrib- 
ution and reproduction in other forums, 
provided the original authors and source 
are credited and subject to any copy- 
right notices concerning any third-party 
graphics etc. 



Frontiers in Physiology | Computational Physiology and Medicine 



February 2013 | Volume 4 | Article 4 | 10 



